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[ A Cylinder Specific Performance 
Parameter Computed for an 
Internal Combustion Engine ] 

Background of Invention 

[000 1 ] Field of the Invention 

[0002] The present invention relates generally to a method for estimating indicated mean 
effective pressure on each cylinder of an internal combustion engine for each firing 
attempt. 

[0003] Background of the Invention 
[0004] 

Computational methods to estimate the power produced in each cylinder of an 
engine have been developed, as described in U.S. Patent 6,223,120. The inventors of 
'1 20 have disclosed a technique in which, based on the torsional characteristics of the 
engine and a measure of engine speed sampled at high bandwidth, the power 
produced in each cylinder can be computed. Such information can be used to 
determine a cylinder providing more or less power than the others. The inventors of 
the present invention have recognized disadvantages in this prior approach because 
the technique, disclosed in '120, computes engine power in each cylinder, but not 
instantaneously. Instead, the method of '1 20 requires sampling over a considerable 
period of time to obtain a reliable measurement. Thus, '1 20 appears to provide a 
measure of engine power in the cylinder only under steady speed conditions and 
recurring conditions. The inventors of the present invention have recognized that an 
intermittent condition, such as an occasional misfire, may not be detected using the 
method described in '1 20. The method in '1 20 is suitable for engines operating at 
steady conditions, such as power generating engines or pumping station engines, but 
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is not well suited for automobile engines in which engine speed changes regularly. 
Another disadvantage is that the method in '120 involves Fourier analysis, which 
would tax onboard engine computers as this type of analysis is computationally 
intensive. 

Summary of Invention 

[0005] Disadvantages of prior art methods are overcome by a method for determining a 
cylinder specific performance parameter for a cylinder of a multi-cylinder internal 
combustion engine, in which a performance parameter representative of engine 
output across a first group of cylinders based on an engine operating parameter and a 
cylinder deviation parameter representative of engine output of a single combustion 
event in the cylinder are determined. Based on the cylinder specific performance 
parameter and the cylinder deviation parameter, the cylinder specific performance 

M 8 parameter can be computed. 
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[0006] The engine rotational acceleration deviation is based on the difference between a 
cylinder specific rotational acceleration and a rotational acceleration average, both 
being determined from the rotational speed of the engine. 

[0007] Preferably, the cylinder specific performance parameter is indicated mean effective 



m pressure (IMEP) for a specific cylinder and said group of cylinders comprise all engine 

Ft S 

? Jf cylinders or a cylinder bank. The engine performance parameter is based on one or 

G more of: air flow rate to the engine, ignition timing, air-fuel ratio, fuel pulse width, 

and a quantity of exhaust gas recirculation. 

[0008] An advantage of the present invention is that an estimate of, or other cylinder 

specific performance parameter, can be obtained under transient conditions as well as 
steady-state. Thus, the method disclosed herein is well suited for engines in which 
speed is frequently changing such as the case for engines installed in automobiles. 

[0009] Another advantage is that an estimate of a cylinder specific performance 
parameter, according to an aspect of the present invention, is much less 
computationally intensive than prior art methods. 



[0010] 



An advantage of the present invention is that it relies on sensors that are available 
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on most current engines. Furthermore, the method does not rely on in-cylinder 
pressure transducers which can be costly and difficult to package. 

[001 1] Yet another advantage is that if the engine were equipped with sensors providing 
redundant measures of the performance parameter (in-cylinder pressure transducers, 
as an example), the estimate of the performance parameter, as provided by the 
present invention, can be used to determine anomalous behavior, sensor drift, or 
aging of the other sensors. Such information can be used to update calibratable 
coefficients in the electronic control unit or to determine when a sensor falls out of 
operational limits. 

[001 2] Although prior art methods are able to detect a misfire in a cylinder, an advantage 
of the present invention is that by providing a measure of a cylinder specific 
performance parameter, such as IMEP, the magnitude of a deficit in a lower 
performing cylinder can be evaluated. Thus, for example, a cylinder which is 
deteriorating over time can be detected prior to the cylinder deteriorating to the point 
of complete misfire. 

[001 3] The above advantages, other advantages, and features of the present invention 
will be readily apparent from the following detailed description of the preferred 
embodiments when taken in connection with the accompanying drawings. 

Brief Description of Drawings 

[001 4] The advantages described herein will be more fully understood by reading an 

example of an embodiment in which the invention is used to advantage, referred to 
herein as the Detailed Description, with reference to the drawings wherein: 

[001 5] Figure 1 is a schematic of an engine equipped with sensors according to an aspect 
of the present invention; 

[001 6] Figure 2 is a graph of measured IMEP for approximately 6000 combustion cycles 
of a four-cylinder engine, measure IMEP is based on in-cylinder pressure transducer 
measurements; 

[001 7] Figure 3 is a graph of AIR_CHG over the same test cycle shown in Figure 2; 
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[001 8] Figure 4 is a graph of engine acceleration computed to determine the contribution 
by individual cylinders for each firing event, taken over the same test cycle as Figure 
2; 

[001 9] Figure 5 is a graph of deviation in engine acceleration for individual cylinders for 
each firing event, taken over the same test cycle as Figure 2; and 

[0020] Figure 6 is a graph of estimated IMEP, according to an aspect of the present 
invention. 

Detailed Description 

[0021] A4-cylinder internal combustion engine 10 is shown, byway of example, in 

Figure 1. Engine 10 is supplied air through intake manifold 12 and discharges spent 
gases through exhaust manifold 14. An intake duct upstream of the intake manifold 
1 2 contains a throttle valve 32 which, when actuated, controls the amount of airflow 
to engine 10. Sensors 34 and 36 installed in intake manifold 12 measure air 
temperature and mass airflow (MAF), respectively. Sensor 31 , located in intake 
manifold 14 downstream of throttle valve 32, is a manifold absolute pressure (MAP) 
sensor. A partially closed throttle valve 32 causes a pressure depression in intake 
manifold 1 2. When a pressure depression exists in intake manifold 1 2, exhaust gases 
are caused to flow through exhaust gas recirculation (ECR) duct 1 9, which connects 
exhaust manifold 1 4 to intake manifold 1 2. Within ECR duct 1 9 is EGR valve 1 8, which 
is actuated to control EGR flow. Fuel is supplied to engine 10 by fuel injectors 26. Each 
cylinder 1 6 of engine 1 0 contains a spark plug 26. A pressure transducer 30 is shown 
installed in each cylinder 1 6. The crankshaft (not shown) of engine 1 0 is coupled to a 
toothed wheel 20. Sensor 22, placed proximately to toothed wheel 20, detects engine 
10 rotation. Engine output shaft 23 coupled to the crankshaft of engine 10 is coupled 
to a transmission (not shown). An in-line torque sensor 35 is coupled to output shaft 
23. Sensor 24, in exhaust manifold 14, is an exhaust gas component sensor. Exhaust 
gas component sensor 24 is an exhaust gas oxygen sensor. Alternatively, exhaust gas 
component sensor 24 is a wide-range oxygen sensor, a nitrogen oxide sensor, a 
hydrocarbon sensor, or other gas component sensor as may become available. 

[0022] 

Continuing to refer to Figure 1 , electronic control unit (ECU) 40 is provided to 
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control engine 10. ECU 40 has a microprocessor 46, called a central processing unit 
(CPU), in communication with memory management unit (MMU) 48. MMU 48 controls 
the movement of data among the various computer readable storage media and 
communicates data to and from CPU 46. The computer readable storage media 
preferably include volatile and nonvolatile storage in read-only memory (ROM) 50, 
random-access memory (RAM) 54, and keep-alive memory (KAM) 52, for example. 
KAM 52 may be used to store various operating variables while CPU 46 is powered 
down. The computer-readable storage media may be implemented using any of a 
number of known memory devices such as PROMs (programmable read-only memory), 
EPROMs (electrically PROM), EEPROMs (electrically erasable PROM), flash memory, or 
any other electric, magnetic, optical, or combination memory devices capable of 
storing data, some of which represent executable instructions, used by CPU 46 in 
controlling the engine or vehicle into which the engine is mounted. The computer- 
readable storage media may also include floppy disks, CD-ROMs, hard disks, and the 
like. CPU 46 communicates with various sensors and actuators via an input/output 
(I/O) interface 44. Examples of items that are actuated under control by CPU 46, 
through I/O interface 44, are fuel injection timing, fuel injection rate, fuel injection 
duration, throttle valve 32 position, spark plug 26 timing, EGR valve 1 8. Various other 
sensors 42 and specific sensors (engine speed sensor 22, in-line torque sensor 25, 
cylinder pressure transducer sensor 30, engine coolant sensor 38, manifold absolute 
pressure sensor 31 , exhaust gas component sensor 24, air temperature sensor 34, 
and mass airflow sensor 36) communicate input through I/O interface 44 and may 
indicate engine rotational speed, vehicle speed, coolant temperature, manifold 
pressure, pedal position, cylinder pressure, throttle valve position, air temperature, 
exhaust temperature, exhaust stoichiometry, exhaust component concentration, and 
air flow. Some ECU 40 architectures do not contain MMU 48. If no MMU 48 is 
employed, CPU 46 manages data and connects directly to ROM 50, RAM 54, and KAM 
52. Of course, the present invention could utilize more than one CPU 46 to provide 
engine control and ECU 60 may contain multiple ROM 50, RAM 54, and KAM 52 
coupled to MMU 48 or CPU 46 depending upon the particular application. 

[0023] 

According to the present invention, a value of a cylinder specific performance 
parameter (estPP.n) for the nth cylinder is computed as the sum of cylinder specific 
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information, DACCEL.n, and information which applies to the engine as a whole, PPall: 
[0024] estPP,n = b * DACCEL.n + PPall (1) 

[0025] where PPall is an engine performance parameter which is determined based on 

output of sensors and/or data available onboard conventional, production vehicles 

and is discussed further below. Constant b is determined as discussed below. 

0 

DACCEL,n is the deviation in engine rotational acceleration of the nth cylinder from 
average engine rotational acceleration (ACCELavg) and, by way of example, is 
computed as: 

[0026] DACCEL.n = ACCEL, n ACCELavg (2) 

where ACCEL, n is engine rotational acceleration of the nth cylinder, which in one 
example is computed according to the method taught in U.S. Patent 5,056,360, 
assigned to the assignee of the present invention and incorporated in its entirety 
herein by reference. In '360, the interval and phase at which engine rotational speed is 
measured is determined to provide accurate and substantially independent rotational 
acceleration measures for each expansion stroke in each cylinder. In a four-stroke 
engine, each cylinder undergoes the processes of intake, compression, expansion 
(combustion), and exhaust to complete a cycle. It is the expansion stroke which 
provides an indication of the efficacy of the combustion process. If the engine is a 
four-cylinder engine, the expansion strokes of the cylinders are evenly spaced over 
two revolutions of engine 1 0. One expansion event occurs every 1 80 degrees and no m 
two expansion strokes overlap in time. It is described in '360 how to determine the 
phase and duration over which to compute ACCEL,n for situations in which the 
expansion stroke of one cylinder overlaps the expansion stroke of another cylinder, as 
is the case in engines with greater than four cylinders. ACCELavg is computed, for 
example, as described in U.S. Patent 5,044,194, assigned to the assignee of the 
present invention and incorporated in its entirety herein by reference. In '194 a sliding 
median filter is used to determine ACCELavg, i.e., ACCELavg is the median of the prior 
measures of ACCEL.n. Alternatively, other measures of average, such as mean, are 
used in place of median. Also, the window over which the prior events are used in 
averaging is determined so as to provide a smooth measure of ACCELavg without 
damping the information about the change in ACCELavg of the engine. The averaging 



□ 



[0027] 



APP_ID=10064192 



Page 6 of 24 



window to determine ACCELavg is on the order of one engine revolution, depending 
on the data bandwidth of crankshaft position from which ACCELavg is computed. 



[0028] The engine performance parameter (PPall) is determined over all the cylinders, 
collectively. It is based on any combination of air charge (AIR.CHG), engine speed 
(RPM), fuel pulse width (FPW), commanded air-fuel ratio (AFR), measured air-fuel 
ratio, concentration of an exhaust gas component, exhaust gas temperature, spark 
advance (SA), intake manifold absolute pressure (MAP), position of an exhaust gas 
recirculation (EGR) valve, accelerator pedal position, barometric pressure, engine 
coolant temperature, fuel injection timing, fuel injection supply pressure, valve timing, 
exhaust backpressure, and others. AIR_CHG is an estimate of the mass of air expected 
to be trapped in the cylinder as a result of an intake stroke. AIR_CHC is computed 
based on the output of the MAF sensor, which provides a measure of the amount of 
air inducted into the engine. Alternatively, AIR.CHC is based on engine RPM and MAP. 
FPW is the duration that a fuel injector is commanded to remain open. Thus, FPW 
relates to the amount of fuel delivered during a fuel injection event. EGR relates to a 
quantity of exhaust gases which are routed to the intake manifold. SA relates to the 
time that the spark plugs are commanded to fire in relation to the position of the 
position. Valve timing relates to the opening and closing events of intake and exhaust 
valves coupled to the cylinders of the engine and refers to engines equipped with 
variable valve timing devices. 

[0029] An example of a cylinder specific performance parameter, which may be 

according to the present invention, is indicated mean effective pressure, IMEP. IMEP,n 
is the amount of work done per cylinder in each cylinder divided by the displacement 
volume of the cylinder: 



where the integral is taken over a complete engine cycle (i.e., 720 degrees 
revolution for a four-stroke engine and 360 degrees for a two-stroke engine), P,n is 
the pressure in cylinder n, V,n is the volume above the piston in the nth cylinder, and 
dV,n is the change in volume. The cylinder pressure, P,n, is measured by cylinder 
pressure transducer 30 in the nth cylinder and V,n may be computed based on engine 
geometry and a measure of the crank angle position of the engine from an toothed 



IMEP,n 




P,n dV,n / DISP 



(3) 



[0030] 
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wheel 20 and sensor 22 coupled to engine 10 or other device, such as an optical 
encoder, for determining crank position. The determination of IMEP, according to 
equation (3) above, is based on pressure measurements from in-cylinder pressure 
transducers 30, typically in an experimental setting. 

[0031] According to the present invention, an estimate of IMEP can be determined 

without relying on in-cylinder pressure transducers 30 by recasting equation (1) in 
terms of IMEP in place of the generalized performance parameter: 

[0032] estlMEP.n = c * DACCEL.n + modellMEPall. (4) 

[0033] An example of a modeled performance parameter, according to the present 
invention is: 

[0034] IMEPall = c i + c * AIR_CHG. (5) 



[0035] Or, a more complicated example: 

- [0036] IMEPall = c ^ + c * AIR.CHG+ c * AIR_CHG + c n * (SA SAMBT)+ 

«♦= 3 4 5 6 

£ [0037] c * MAP + c * SA * AIR.CHG + c / AFR + c * AIR.CHG 2 . (6) 
iU 7 8 9 10 



[0038] quation 6 includes three engine parameters with linear, squared, cross terms, etc. 
SAMBT is MBT spark advance, which is the spark advance which would provide the 
q highest torque. SAMBT is a known quantity for a given set of conditions. Alternatively, 

U modellMEPavg is determined via a neural network model, look-up tables, a physics 

based model from first principles, or other models. 

[0039] By substituting equation (4) into equation (5), e 

[0040] estlMEP.n = c * DACCEL,n + c 1 + c ^ * AIR.CHG. (7) 



[0041] 



To determine the values of the constants, c , c i , and c ^ , data are collected via 
laboratory grade pressure sensors installed in each cylinder of engine 10. Cylinder 
specific pressure data are collected at the same time that other engine parameters are 
collected. In the case of equation (7), MAF is measured and recorded, from with 
AIR_CHG can be computed. Alternatively, if the model were based on equation (6), 
MAP, MAF, SA, and AFR would be measured and recorded. The data are regressed, 
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inserting the measured IMEP in equation (7) as the term estlMEP,n, to determine the 
values of the constants. 

[0042] Additionally, the measured IMEP data via the laboratory pressure measurements 
are used to validate the form of the model employed, i.e., the appropriate form of 
estlMEP,n. That is, to determine whether a simple form of the model, such as equation 
(6), or a complicated form of the model, such as equation (7), is more appropriate. 

[0043] Finally, if the vehicle were equipped with in-cylinder pressure transducers, the 

present invention can be used to advantage to detect errors in the measurement due 
to sensor drift, sensor malfunction, electrical integrity, or others. 

To determine the efficacy of the present invention, data was collected from a four- 
cylinder, homogeneous-charge, spark-ignition engine, which was operated at a 
variety of operating conditions and with intentional misfiring in one cylinder. 
Laboratory grade in-cylinder pressure transducers 30 were used to collect cylinder 
pressure data from all four cylinders. IMEP was determined based on pressure 
transducer measurements in each cylinder. In Figure 2, the resulting IMEP is shown for 
each sequential firing for cylinders 1 , 3, 4, and 2 (firing order of this particular engine) 
as they undergo approximately 1 500 combustion events each (6000 total for the 
engine). The engine is being operated at a fairly low IMEP condition between events 0 
and 900. From events 900 to about 1 500, IMEP is at a moderate level. At 

th 

approximately event 2000, a misfire was caused to occur in approximately every 25 
event by interrupting spark plug firing every 25 events. Consequently, every time that 

a spark plug fails to fire, that cylinder has a negative IMEP indicating that the misfiring 

cylinder absorbs more energy in undergoing a cycle than it provides. The data 

represented in Figure 2 are derived from laboratory grade in-cylinder pressure 

transducers 30 and laboratory quality data acquisition systems. Having such data 

available on board a production vehicle is highly desirable. However, production 

engines are rarely equipped with in-cylinder pressure measuring equipment. A 

primary aspect of the present invention is to provide an alternative to in-cylinder 

pressure measuring equipment to obtain substantially similar information as that 

represented in Figure 2. 

[0045] As men tioned above, the present invention relies on production sensors or data 



O 



fli 

6 



[0044] 
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available within ECU 60. Data based on one such sensor, a MAF sensor, is shown in 
Figure 3. The data of Figure 3 are over the same interval as those shown in Figure 2. 
Throttle valve 32 was adjusted during the test period causing AIR_CHG of the engine 
to vary by a factor of three over the test period. By comparing Figure 3 to Figure 2, it 
can be seen that the general shape of the two are in excellent agreement. However, 
because AIR_CHG is a quantity that applies to the entire engine, i.e., no cylinder 
specific information is contained therein, some of the features of Figure 2 are missing 
in Figure 3. Most notably, the misfires that are shown in Figure 2 between events 
2000 and 6000 are not detected in Figure 3. The present example utilizes only 
AIR_CHC information as shown in equation (5). 

[0046] Engine rotational acceleration (ACCEL) is plotted in Figure 4, shown for the same 

test duration of Figures 2 and 3. The rotational acceleration is computed in such a way 
u. as to correspond with the combustion events of each of the cylinders. Figure 4 shows 

D consecutive measurements of the contribution of cylinders 1,3,4, and 2 in 

O 

fTi succession. ACCEL, shown in Figure 4, is computed based on a signal from a 

production engine speed sensor, such as 22 shown in Figure 1 . DACCEL is computed 
=0 based on ACCEL and shown in Figure 5. 

5=5 = 

\ y 

L [0047] In Figure 6, estlMEP,n is shown, which is computed according to equation (7) 

u * 

01 above: 

ru 

g [0048] estlMEP.n = c * DACCEL,n U ] + c ^ * AIR_CHG (7) 

n i 

[0049] where c q( c ] , and c ^ are found by regressing measured IMEP data, computed, 
according to equation (7) above, against DACCEL.n and AIR.CHG, where the measured 
IMEP is inserted in for estlMEP,n. After determining the values of the constants, 
estlMEP.n was computed based on production sensors (a MAF sensor in the present 
example) and is shown in Figure 6. The results in Figure 2 based on laboratory grade 
equipment and the results in Figure 6, according to the present invention, are in 
excellent agreement. 



[0050] 



The efficacy of the method has been shown for the test conditions employed 
represented in Figures 2-6. Figure 6, which was determined according to an aspect of 
the present invention, excellently represents the data collected by in-cylinder pressure 
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transducers, which is represented in Figure 2. The success depends, in part, on the 
choice of the model employed, which in the present example is equation (5). To 
satisfy unusual operating conditions in a homogeneous-charge, spark-ignition 
engine, a more complicated form of the engine performance parameter (PPall) may be 
found to be more accurate than equation (5) used for the example situation discussed. 
Additionally, the present invention can be applied to engines with other combustion 
types, such as: ldiesel, homogeneous-charge compression ignition, direct injection 
spark ignition, and others. Alternatively, models based on other data, eg., FPW, MAP, 
ECR, AFR, and combinations thereof, yield a more suitable, or accurate, model of PPall. 
As mentioned above, other types of models, such as neural nets, lookup tables, first 
principle physics based models, are other examples of models that can be used to 
describe any of the engine types discussed. 

[0051] Alternatively, instead of using D ACCEL, n in the above equations, which is the 
5 deviation in engine rotational acceleration of the nth cylinder from average engine 

gi rotational acceleration, a torque deviation of the nth cylinder, DT,n, is used. The 

signal from in-line torque sensor 25 is used to determine torque deviation associated 
O with the torque delivered by each of the cylinders. It is known by those skilled in the 

s art that it is difficult to obtain a reliable measure of absolute torque from an in-line 

y torque sensor. However, in-line torque sensors accurately measure deviations in 

m 

nj torque, which is suitable for use in the present invention. Thus, to use DT,n in place of 

Q 

S DACCEL,n, a suitable model is: e 



[0052] stlMEP,n = c ' * DT,n + C ] + c * AIR.CHG (8) 

[0053] where c Q ' is a calibration constant determined analogously to c of equation 7. 
Alternatively, other models involving DT,n may be found to more accurately represent 
estlMEP.n. Furthermore, another alternative is to base estlMEP,n on both DT,n and 
D ACCEL, n if both an inline torque sensor and a speed sensor are available. 

[0054] The specific example provided has estlMEP.n as the performance parameter. 
Alternatively, work done per cylinder or torque contribution per cylinder can be 
computed analogously as IMEP. The above equations may be recast for these 
alternative quantities or other quantities based on combinations of IMEP, work, and 
torque. 
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[0055] Engine performance parameter has been described in terms of a quantity, such as 
indicated mean effective pressure, that applies to the engine as a whole. However, 
there are situations in which it is useful to define a performance parameter which is 
associated with a group of cylinders, where the group of cylinders is less than all 
cylinders in the engine. For example, a V-8 engine comprises 2 four-cylinder banks of 
cylinders. There are occasions in the operation of the engine that it is desirable to 
command a different spark advance, fuel-air ratio, or other parameter to each bank of 
the cylinder. In this case, one or more performance parameter can be determined for 
one or more groups of cylinders. Another example is in an engine with a variable 
displacement engine (VDE). The displacement is varied in a VDE by deactivating one or 
more cylinders in the engine, thus operating on a lesser displacement than the 
engine's displacement as a whole. Typically, valves are deactivated to deactivate a 
cylinder. In this case, a performance parameter is computed for the group of cylinders 
that are not deactivated, i.e., activated cylinders. 

[0056] While several modes for carrying out the invention have been described in detail, 
those familiar with the art to which this invention relates will recognize alternative 
designs and embodiments for practicing the invention. The above-described 
embodiments are intended to be illustrative of the invention, which may be modified 
within the scope of the following claims. 
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